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ABSTRACT 

We have used the Green Bank Telescope to measure H I absorption against 
the anomalous X-ray pulsar XTE J1810-197. Assuming a flat rotation curve, we 
find that XTE J1810-197 is located at a distance of 3.4^Q7kpc. For a rotation 
curve that incorporates a model of the Galactic bar, we obtain a distance of 
4.0lo'8 kpc- Using a rotation curve that incorporates a model of the Galactic bar 
and the spiral arms of the Galaxy, the distance is 3.7 ± 0.6 kpc. These values 
are consistent with the distance to XTE J1810-197 of about 3.3 kpc derived from 
its dispersion measure, and estimates of 2-5 kpc obtained from fits to its X-ray 
spectra. Overall, we determine that XTE J1810-197 is located at a distance of 
3.5 ±0.5 kpc, possibly not far in front of the infrared dark cloud GlO.74-0.13. We 
also used the GBT in an attempt to measure absorption in the OH ^vr3/2(J = 3/2) 
lines against XTE J1810-197. We were unsuccessful in this, mainly because of 
its declining radio fiux density. Analysis of H I 21cm, OH ^7r3/2(J = 3/2), and 
^^C0(2 — s> 1) emission toward XTE J1810-197 allows us to place a lower limit of 
A^H ^ 4.6 X 10^^ cm~^ on the non- ionized hydrogen column density to XTE J1810- 
197, consistent with estimates obtained from fits to its X-ray spectra. 

Subject headings: ISM: clouds — ISM: individual (GlO.74-0.13) — pulsars: in- 
dividual (XTE J1810-197) — radio lines: ISM 



Introduction 



Prior to the de tection of puls e d rad io emission from the anomalous X-ray pulsar (AXP) 
XTE J1810-197 by ICamilo et al.l (120061 ) . pulsed emission from the dozen known magnetars 
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had been detected in X-rays in all instances, and in one case at optical wavelengths. Es- 
timating the distance to a magnetar has relied on associating it with a supernova remnant 
(SNR) of known distance, or by fitting its X-ray spectrum and parameterizing the energy- 
dependent absorptio n by interstellar gas along the line of sight with a non-ionized hydrogen 
column density iVn (IMorrison fc McCammonlll983l ). Using standard relations, A^h is related 
to visual extinction Ay, whic h is turned into a di stance estimate using mean values in the 
Galactic plane (see, e.g., §6 of lGotthelf et al.ll2004l ). o r is directly calibrated as a fun ction of 
distance using stars of known luminosity in the field (IDurant fc van Kerkwijkl 120061 ). Some- 
times, the probab le location of the X -ray source in a well-studied star cluster can be used to 
infer its distance (IMuno et al.ll2006l ). 



The unique detection of pulsed radio emission from XTE J1810-197 allows one to deter- 
mine its distance using methods that are not applicable to other magnetars. The dispersion 
measure (DM, the total column density of free electrons alon g the line of sight) was obtained 
upon discovery of the radio pulsations (jCamilo et al.l 120061 ). A model for the Galactic free 
electron distri bution then yields a di stance estimate, in this case d ^ 3.3 kpc using the most 
recent model (ICordes fc Lazidl2002l ). Also, bright, pulsed radio emission allows kinematic 
distance limits to be obtained by observing spectral lines that are seen in absorption against 
the magnetar, as we report here using H I. 

Obtaining a reliable distance to XTE J1810-197 allows for a precise determination of the 
luminosity of the star based on its measured flux in a variety of wavebands. The distance also 
allows a proper motion (see lHelfand et al.l 120071 ) to be converted into a tangential velocity. 
Magnetars are thought to be very young neutron stars, and are expected to be found near 
star forming regions and/or spiral arms. Knowing the distance to XTE J1810-197 allows 
this prediction to be tested in this case. Besides providing a kinematic distance, the H I 
absorption spectrum can also give an independent estimate of N^, which may be compared 
to results from X-ray spectral fitting. 

In § 2 we present the H I and OH observations and data analysis. This is followed by 
a determination of the hydrogen absorption spectra toward XTE J1810-197, in § 3, and of 
its kinematic distance in § 4. In § 5 we comment briefly on some features of the neutral 
hydrogen toward XTE J1810-197, and in § 6 on the OH absorption limits. We obtain a limit 
on the hydrogen column density in § 7, and comment on models of the line of sight toward 
XTE J1810-197 in 5 8. We conclude in 5 9 with a discussion of our main results. 
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2. Observations and Data Analysis 

2.1. HI 21cm Absorption Observations 

XTE J18 10-197 was observed with the National Radio Astronomy Observatory (NRAOQ) 
Robert C. Byrd Green Bank Telescope (GET) for approximately 3 hr on 2006 June 6 in order 
to measure the H I 21 cm line absorption. The GET has an unblocked aperture, a spatial 
resolution of 9'.2 at 21 cm, and a system temperature on cold sky of 18 K. The NRAO spec- 
tral processor (SP), an FFT spectrometer, was used as the detector. The SP provides good 
dynamic range for the observations via its 32-bit sampling, and was used in its pulsar mode 
whereby it accumulates spectra that are folded synchronously with the pulsar period. The 
SP was configured to have 1024 channels for each linear polarization across a bandwidth 
of 2.5 MHz, producing a spectral resolution of 0.52 km s~^ per channel. We recorded 128 
spectra (phase bins) evenly spaced in time across each individual period of XTE J1810-197. 
The duty cycle of this 5.54 s pulsar was such that 4-5 of these bins contained pulsed flux. 

2.2. OH ^7r3/2(J = 3/2) Absorption Observations 

XTE J1810-197 was observed with the GET for approximately 38 hr in order to measure 
absorption in the 1612, 1665, 1667 and 1720 MHz ^7r3/2(J = 3/2) ground state transitions of 
OH. Observing details are given in Table [H The GET's angular resolution is approximately 
8' for each of the OH transitions. The spectral processor was also used for these measure- 
ments and was conflgured to provide 256 spectral channels for each linear polarization across 
each of the four OH transitions. A bandwidth of 0.625 MHz was used, giving a spectral reso- 
lution of 0.44 km s~^ per channel. We chose this narrower bandwidth (with higher frequency 
resolution) since we already knew which velocities displayed H I absorption. As with the H I 
measurements, the pulsed flux of XTE J1810-197 was detected in 4-5 phase bins. 

The GET auto-correlation spectrometer (ACS) was used to obtain pulsar "off" spectra 
for the four OH lines. The data were obtained in six position-switching observations, each 
consisting of two minutes on source followed by two minutes off source, resulting in an 
effective integration time of 12 minutes. An off position two minutes of time in R.A. offset 
from the position of XTE J1810-197 was used so that approximately the same hour-angle 
coverage was obtained for the on and off positions. The ACS was conflgured to observe all 
four of the OH ^7r3/2(J = 3/2) lines simultaneously with 8192 channels within a bandwidth 
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of 12.5 MHz, resulting in a spectral resolution of 0.27kms ^. These observations allow for a 
better calibration of the pulsar "off" OH spectra. 



2.3. Data Reduction 



The data we r e ana lyzed using a method similar to that described in IWeisberg] (Il978l ) 
and lMinter et al.l (120051 ). For each individual pulse of XTE J1810-197 and for each spectral 
line we measured the flux T(z/, 0, V) in terms of the equivalent brightness temperature T, 
versus frequency u, phase across the pulse (equivalent to time), and polarization V. For 
each polarization and pulse, the data were divided into two separate parts, a pulsar "on" 
spectrum and a pulsar "off" spectrum: 



(T(^,<Aon,P)).(„o,bs) 



(^(^.'^on.^)).(noabs) 



(1) 



N, 



(2) 



4'oft 



where (pon are the phase bins when the pulsar is on, (pos are the phase bins when the pulsar 



is off, N(f,^g is the total number of phase bins when the pulsar is off and ()j,(i™bs) denotes 
averaging over the frequencies that do not show absorption in the final spectruno. A pseudo- 
absorption spectrum for the ith pulse is then formed by taking the difference between the 
pulsar on and pulsar off spectra. 



T^,{u,V)-T:^{u,V) = {T;{u,V)e-^^''^+Tu iiu)+T,ys) - (Th i(z/) +T, 

= T^(i/,P)e-"(") 



sysy 



(3) 



where Tp is the brightness temperature of the pulsar. We then take a weighted average of 
the pseudo-absorption spectrum 



(Tj5(!/,-p))_^(noabs) 



(4) 



^An iterative approach in the data reduction is necessary in order to determine z^(noabs). 
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and store the weights for later use when the polarizations are averaged together. This 
weighting is proportional to the signal-to-noise ratio obtained for each pulse of XTE J1810- 
197. A third-order orthogonal polynomial was fitted to the {Tp{u, V)) e"'^'-^^ spectrum in order 
to determine the intrinsic pulsar brightness, Tp*(z/, "P), at all frequencies, i.e., by extrapolating 
(Tp(i/, P)) across the absorption features. This yielded the absorption spectrum 

for each polarization. The absorption spectra for the different polarizations were then aver- 
aged using the weights from equation (jlj) to create the final, measured absorption spectrum 
for each spectral line. 

For the H I data, the pulsar off spectrum (i.e., the normal H I emission spectrum) 
was converted from detector counts to a Kelvin scale using a calibrated noise diode that 
was injected during two-minute calibration scans. Observations of the lAU H I standard 
source S6 were also made to put the pulsar off spectrum on the lAU standard brightness 
temperature scale. 



3. The Absorption Spectra Toward XTE J1810-197 

The H I absorption spectru m toward XTE J1810 -197 is shown in Figuredl XTE J1810- 



197 is highly linearly polarized (jCamilo et al.ll2007bl ) so that the YY polarization signal had 



a much better signal-to-noise ratio (by a factor of ~ 4) than the XX polarization. The two 
polarizations were averaged together with weights given by their signal-to-noise ratios to 
produce the spectrum shown. The spectrum has not been smoothed in any way and shows 
the native resolution of the observations. 

Five Gaussians were fitted to the opacities determined from the H I absorption toward 
XTE J1810-197. The results of the fit are shown in Table [2] and in Figure [2j The columns in 
Table [2] give for each line, respectively, the opacity. Local Standard of Rest velocity (^sr); 
and full- width at half-maximum (FWHM). The errors listed in Table [2] are la errors from 
the Gaussian fits. 



The Vlsr = 7.7kms ^ absorption line is associated with the Hee schen Cloud (jRiegel fc Crutcher 



19721) . which has also been called the Riegel-Crutcher Cloud by iMcClure-Griffiths et al. 
km^ . The Heeschen Cloud is a nearby, d ^ I2h ± 25 pc, cold cloud (T^pin ^ 40 K) that 
covers Galactic longitudes 345° to 25° and latitudes ±6°. This cloud is also seen as a self- 
absorption feature in the H I emission spectrum (top panel of Fig. [1]). The weak absorption 
line at Vlsr = 14.1 kms~^ can be kinematically associated with gas in the Carina-Sagittarius 
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spiral arm assuming that th e line arises on the nea r side of the tangent point and using the 
Galactic rotation model of lEnglmaier fc GerhardI (120061 ) . Likewise, the Vlsr = 22.8 and 
25.7 km s^^ lines can be kinematically associated with the Crux-Scutum spiral arm. 

The OH absorption spectra toward XTE J1810-197 are shown in Figures [3H6l With 
the strength of the pulsar emission during the OH measurement epochs being much weaker 



than was the case for the H I measurement epoch (jCamilo et al. 



2007al ). the strong linear 



polarization of XTE J1810-197 means that only one polarization effectively contributed to 
the measurement of the OH absorption against XTE J1810-197. Thus, the OH absorption 
data shown in Figures [3H6] and discussed in this paper are only from a single linear polariza- 
tion (YY). The la opacity limits for any OH absorption toward XTE J1810-197 are given 
in Table El 



4. Determining the Kinematic Distance to XTE J1810— 197 



4.1. Near or Far Side of the Velocity— Distance Relationship? 

XTE J1810-197 is located at Galactic coordinates {l,b) = 10?726, -0?158. At this 
Galactic longitude the kinematic velocity-distance relationship is double- valued (see Fig. [7]). 
Our first concern is whether it is possible to determine on which side of the tangent point 
XTE J1810-197 lies. At this Galactic longitude, the tangent point is abou t 8.3 kp c from the 
Sun and has a ^sr ~ 168kms~^ using the flat rotation curve of iBurtonI ( 19881). Distanc e 



estimates for XTE J1810-197 based on its X-ray-fitted are 2.5-5 kpc (iGotthelf et al. 



2004 : iGotthelf fc Halpernll20(35l : IPurant fc van KerkwiikI 120061 1 . and the DM-based distance 
is 3.3 kpc (ICamilo et al.l 120061 ). These strongly suggest that XTE J1810-197 lies on the near 
side of the tangent point. 

Another magnetar, SGR 1806-20, which lies 41' from XTE J1810-197, has H I absorp- 
tion at velociti es greater than Vt.rr = 50kms~^ ari d is thought to lie on the far side of the 
tangent point jMcClure-Griffiths fc Gaensleil l2006h. Both SGR 1806-20 and the Galactic 
SNR Gil. 2-0. 3 (IBecker et al.l 119851 ) also have H I absorption features at negative velocities 
that must arise from the far side of the Galactic disk if the rotation curve is fiat. However, 



as the Galactic rotation model of IWeiner fc SellwoodI (119991 ) shows, these negative velocities 
can also arise in the Galactic bar at a distance of 7 kpc toward these sources. The Galactic 
bar induces large radial motions that deviate from the normally assumed circular Galactic 
rotation. Along the line of sight to XTE J1810-197, the Galactic bar is responsible for 
motions > 90kms~^ more negative than the prediction of a fiat rotation curve (see Fig. 
Negative-velocity H I absorption is not found toward XTE J1810-197, which according to 
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this model establishes an upper-limit to its distance of 7kpc (i.e., it is closer than the Galac- 
tic bar). This strengthens the conclusion that XTE J1810-197 lies on the near side of the 
tangent point. 



4.2. Is the Last Absorption Feature at the Actual Distance or a Lower Limit? 



The number of H I absorption features per kpc in the inner Galaxy is about 1-3 
([Garwood fc Dickeyl Il989l . see their Table 3). Assuming a flat rotation curve and a dis- 
tance of 5 kpc, corresponding to a maximum Vlsr ~ 40kms~^, along with a line width of 
3kms~^ for the average H I absorption feature, we expect that 38%-100% of the velocities 
within Vlsr = 0-40 km s~^ should contain H I absorption along the line of sight toward 
XTE J1810-197. If ~ 38% of the velocity space were occupied by absorption features, then 
the highest velocity feature in the absorption spectrum toward XTE J18 10-197 should be 
taken as indicating a lower limit for the distance, since there is likely a significant distance 
between XTE J1810-197 and the H I-absorbing cloud nearest it. If on the other hand the 
whole velocity range should contain H I absorption, it is likely that the last absorbing feature 
gives the actual distance to XTE J1810-197. 

It is useful to look at absorption measurements toward objects near XTE J1810-197 
in order to help us determine the H I absorption velocity coverage. The line of sight to 
XTE J1810-197 lies near that to the giant Galactic H II complex W 31 (see Fig. 1 of 
Corbel fc Eikenberryl |2004| ). W 31 is comprised of several individual H II regions includ- 
ing GlO.2-0.3, GlO.3-0.1, and GlO.6-0.4. Roughly on the opposite side of W 31 from 
XTE J1810-197 lies GlO.0-0.3, the wind-blown bubble of LBV 1806-20, and SGR 1806-20. 
At higher Galactic longitude than XTE J1810-197 lies the SNR Gil. 2-0.3. All of these 
objects have H I absorption measure ments, except for GlO.0-0. 3 and LBV 1806-20, which 
have NH3 absorption measurements (ICorbel &: Eikenberryl |2004| ). 



Becker et al.l (119851) claimed th at SNR Gil. 2-0. 3 lies well on the far side of the tangent 
point. However, I Green et al.l (Il988l ) noted the lack of absorption between Vlsr = 45kms~^ 
and the tangent point velocity, which led them to assign a distance of ~ 5 kpc to G11.2- 
0.3, while attributing weak absorption at negative velocities to peculiar motions in local 
gas. A dist ance of ^ 5 kpc is also in reasonable agreement with the expected expansion 



of the SNR (IGreen et al. 



19881 ). If we take the Galactic bar into account using the model 



of IWeiner &: SellwoodI (11999), we can reconcile the H I absorption with the expansion/age- 



derived distance estimate of iGreen et al.l (119881 ) . The H I absorption to Gil. 2-0. 3 at I^sr ~ 
— 20kms~^ arises from the Galactic bar, which is at a distance of ~ 5.5-7.0 kpc along this 
line of sight (see Fig. [S]). With both the Galactic bar and Gil. 2-0. 3 being on the near side 



- 8 - 



of the tangent point toward Gil. 2-0. 3, H I absorption is not expected at Vlsr > 45kms ^. 
So we judge that Gil. 2-0. 3 hkely hes within the Galactic bar at ^ 5.5-7kpc. 



Corbel fc Eikenberrvl (120041 ) determined that GlO.3-0.1, GlO.0-0.3, LBV 1806-20, and 



SGR 1806-20 all lie on the far side of the tangent point, while GlO.2-0.3 and GlO.6-0.4 lie on 
the near side. The H I absorption o f GlO.3-0.1 (Fig. 2 of lKalberla et al.lll982h and SGR 1806- 
20 (Fig. 2 of ICameron et al.lboO^'ih fill the range Vlsr = 0-40 km s'^ H I toward GlO.6-0 4 
shows absorption at all velocities within Vlsr = 0-45kms~^ fFig. 32 oflCaswell et al.lll975l) . 
Likewise, the H I absorp tion toward Gl 0.2 -0. 3 (Fig. 4 of Greisen &: Lockman llQTol ) and 



SNR Gil. 2-0. 3 (Fig. 3 of lBecker et al.lll985l ) completely cover the velocity range Vlsr = 0- 
45kms-^ Since GIO. 6-0.4, GlO.2-0.3 and Gll.2-0.3 all lie on the near side of the tangent 
point and are only 16', 33' and 30', respectively, in projection from XTE J1810-197, we 
argue that all velocities between the Sun and XTE J1810-197 should show H I absorp- 
tion. Therefore, the highest-velocity absorption feature seen along the line of sight toward 
XTE J1810-197 gives us its true distance, rather than a lower limit. 



4.3. The Kinematic Distance to XTE J1810-197 



For the kinematic distance models that we discuss we will assume that the Sun is 
located 8.5 kpc from the Galactic center. We will also use a velocity of = 220kms~^ 
as the azimuthal velocity of the LSR about the Galactic center. These are the current 
lAU standard s. Howeve r , there is evidence that the se values may not be correct (see the 
discussions in iReidI 1 19931 : lEnglmaier fc Gerhard! |2006| ) , which would require a scaling of the 
kinematic distance estimates presented below. 

In Figure [7] we show the deter mination of the distance to XTE J1810-197 using the 
flat rotation curve of iBurtonl (119881 ). Cold H I clouds have random motions superposed on 
the uniform Galactic rota tion as indicated by their measured velocity dispersion of 7kms~^ 
(ILockman fc Dickeylll990l ). The dotted lines on either side of the rotation model in Figure [7] 
indicate the deviation from Galactic rotation of ±7kms~^. (The random motions of the 
clouds are relative to the Galactic rotation and not to the measured velocity of the cloud. 
Thus, the 7kms~^ represents an error in converting the model's distance into a velocity and 
not in using the measured velocity of the cloud to obtain a distance as is usually presumed. 
This distinction can be quite important in determining the distance and its errors when the 
observed velocity is at or near a local maximum or minimum in the model. However, in most 
cases it makes only a minor difference.) The arrow at a constant ^sr in Figure [7] indicates 
the velocity of the highest velocity H I absorption. The shaded regions in Figure [7] indicate 
the allowed kinematic distances from the flat rotation curve. Figures [S] and M are similar 
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to Figure [71 using ins t ead t he Galactic rotation models of IWeiner fc SellwoodI (Il999l ) and 
Englmaier fc Gerhard ( 2006 1^. respectively. 



For the flat rotation curve used in Figure [3, the last H I absorption feature toward 
XTE J1810-197 at Vlsr = 25.7 km s~^ (see Table [2]) gives a kinematic distance of 3.4lQ;7kpc 
or IS.Sil^ogkpc. We can rule out the larger distance since this is on the far side of the 
tangent point. Although a fiat rotation model may provide a reasonable distance estimate 
for XTE J1810-197, there are Galactic structures along the line of sight that also need to 
be taken into accou nt. Spiral arna s can have a substantial effect on the expected Galactic 
rotational velocities (lRobertslll972l ). The Galactic bar may also have an influence by adding 
non-circular motions to the general Galactic rotation. 



The rotation model of IWeiner fc SellwoodI (Il999l ) used in Figure E] takes into account 
the effects of a strong potential associated with the Galactic bar. This model uses a global 
Galactic gravitational potential to determine the motions of the gas in the ISM. It uses both 
the minimum and maximum observed H I velocities along lines of sight in the inner Galaxy 
in fitting the properties of the Galactic potential. We can again immediately rule out the 
distance on the far side of the tangent point. Since no H I absorption is seen at velocities 
between 30 and 40kms~^ (see Table [2]) we can eliminate all distances except 4.0lQ'gkpc. The 



Weiner fc SellwoodI Jl999h model, however, does not take into account the effects of spiral 
arms. Since magnetars are expected to lie in or near spiral arms this could have a significant 
impact on the velocity-determined distance to XTE J1810-197. 



The rotation model of lEnglmaier fc GerhardI (120061 ) used in Figure [H takes into account 
the effects of a p otential associated with the Galactic bar as well as those associated with 
spir al arms. The lEnglrnaier &: GerhardI (120061 ) model w as determined in a sirn il ar fas hion as 
the IWeiner fc SellwoodI Jl999h model. However, in the lEnglmaier fc GerhardI (120061 ) model 
only the maximum (minimum) observed velocities for CO were used for positive (negative) 
Galactic longitudes. It should be noted that CO observations do not trace column densities 
as low as H I do es. Combined with fitt i ng to only one side of the velocity-longitude data, this 
means that the lEnglmaier fc Gerhard ( 20061) mod el uses a weaker po tential for the Galactic 
bar co mpared with the lWeiner fc SellwoodI (119991 ) model. In fact, the lEnglmaier fc Gerhard 
( 120061 ) model predicts that there should be no negati v e velo city gas along the line of sight 
toward XTE J1810-197, while the lWeiner &: SellwoodI (Il999l) model does. A s can be seen in 
the Southern Galactic Plane Survey H I data (IMcClure-Griffiths et al.ll2005n. there is plenty 
of gas at negative velocities along this line of sight. We sug gest that the lEnglmaier fc Gerhard 
(120061 ) Galactic bar potential is too weak and that the IWeiner fc SellwoodI (119991 ) model 



■^We have scaled the Englmaier fc Gerhard ( 20061 ) model for Rq = 8.5 kpc 
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should be used for velocities and distances associated with the Galactic bar. 



We can, however, still use the lEnglmaier &: GerhardI (120061 ) model to determine a dis- 
tance for XTE J1810-197 since it i s like ly located in front of the Galactic bar as is ev- 
ide nced bv the IWeiner fc SellwoodI (119991 ) mode l results from abov e . In fact, comparing 
the Englmaier fc Gerhard ( 20061 ) result with the Weiner fc Sellwood (1999) result provides 
an indication of how spiral arms may be affecting the velocity-distance relationship toward 
XTE J1810-197. We can again rule out distances that lie on the far side of the tangent point 
and distances that include 30-40 km s~^ gas along the line of sight. As can be seen from 
Figure [9l we obtain kinematic distances of 3.3]'^Q;2kpc or 4.0l;o;4kpc. Due to uncertainties 
in the modeling of the Galactic rotation curve, the 0.1 kpc gap in distance between these 
two values is in effect negligible. We thus combine these two possible distance ranges into a 
single value, 3.7 ± 0.6 kpc. 

All of the above Galactic rotation models are derived from empirical models fit to 
observed data. Ther e is one model that is fully derived from observational data, that of 
Brand fc Blita (119931 ). In this model, observations of H II regions are used. The velocities 
of the H II regions are determined from recombination lines and H I absorption. Inde- 
pendent measuremen ts provide distances to the H II regions. Using the rotation model of 
Brand &: Blita (119931 ) . shown in Figure [TOl we obtain d = 2.4 ±0.5 kpc. However, this model 
does not have enough H II regions at > 2 kpc in the gen eral direction of XTE J1810-197 
to be able to provide a good velocity-distance relation (see lBrand &: Blitzlll993l ). so we give 
it little weight. 

We conclude that the best estimate of a kinem atic ally determined distance to XT E J1810- 
197 is provided by the lWeiner fc SellwoodI (119991 ) and lEnglmaier fc GerhardI (120061 ) models, 
giving d = 3.1-4.3 kpc. 



5. The Neutral Hydrogen Toward XTE J1810-197: GBT Versus SGPS 



Th e GBT has a resolution of 9. 2 at 21cm while the Southern Galactic Plane Survey 
(SGPS. lMcClure-GrifEths et al.ll2005l ) has a resolution of about 3'.3. In Figure [TT1 we compare 
the GBT H I spectrum with the SGPS spectrum. For the H I self-absorption associated with 
the Heeschen Cloud we find the remarkable result that the line width of the absorption 
increases with increasing spatial resolution! To our knowledge this effect has never been 
observed before for H I on arc-minute resolutions. 



Inspection of the SGPS H I data cube shows that there is structure within the H I 
emission inside of the GBT beam. We convolved the SGPS data with a beam the size of 
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GBT's, so that the two data sets would have the same spatial resolution. As can be seen from 
Figure [TT] the GBT spectrum and 9'.2 resolution SGPS spectrum are in very good agreement. 
If we compare the line width of the H I absorption due to the Heeschen Cloud at ~ 8kms~^ 
with the SGPS data having H I absorption seen against XTE J1810-197, we find that the 
line widths are identical. 

These properties suggest that there is definite spatial structure in the H I emission on size 
scales of 3f3 to 9'.2 (0. 12-0.33 pc for a Heeschan Cloud distance of d = 125 pc). That the tiny 
beam of absorption against the magnetar (limited by the size of the pulsar emission region, 
which is smaller than the pulsar's light- cylinder radius, or about 1 light-second) has the same 
line width as the 0.12pc-wide beam of the SGPS at the Heeschen Cloud, suggests that there 
are few if any spatial structures left unresolved by the SGPS in the Heeschen Cloud. Since 
the H I self-absorption line width does change between the GBT and SGPS resolutions we 
can infer that either the absorption feature comprises multiple narrow line features or that 
the non-thermal broadening of the absorption feature changes between different structures 
in the Heeschen cloud. Attempts to fit the different line-width absorption-line structures in 
the Heeschen cloud with multiple Gaussian components does not improve the fitting, which 
suggests that a single Gaussian is sufficient. This implies that the non-thermal contribution 
to the line width varies within the cloud. 

The standard assumption is that non-thermal line broadening arises from turbulence. If 
the turbulence is intermittent, we can expe ct that the turbulence has decayed more in some 



places than in other locations (lFrischlll995l ). This can then easily explain the observed line 
widths of the Heeschen Cloud H I absorption. The areas that have undergone more damping 
of the turbulence will have less turbulent energy and thus have narrower non-thermal line 
widths. 



6. OH Absorption Limits Toward XTE J1810-197 



Although we did not detect any OH absorption against XTE J1810 -197, the limits are 
still interesting. All previous OH absorp tion detections against pulsars (IStanimirovic et al. 



ff 

20031 : iWeisberg et al.l 120051 : iMinterl 120051 ) have found that the absorption is deeper (larger 
opacity) and has a narrower line width than the pulsar "off" spectra: in the three known 
cases, the OH absorption against the pulsar is 2-3 times deeper than seen in the pulsar off 
spectra. 

For our OH absorption limits against XTE J18 10-197 of r < 0.1 (1 cr), we would expect 
our pulsar off spectra to have OH opacities < 0.033. In Table H] we hst the opacities of the 



- 12 - 



OH lines in the pulsar off spectra. We subtracted the system temperature from the raw 
OH spectra, and then fitted a third order polynomial to determine the continuum emission 
levels. Dividing the spectra by the continuum emission results in e""^ spectra for the pulsar 
off spectra. 

If we assume that the OH pulsar off spectra should have a factor of 2-3 weaker opacity 
than the OH absorption against the pulsar, then we were likely within a factor of about 
2 of detecting OH absorption against XTE J1810-197 at velocities within the range ob- 
served for H I absorption (Table [2]), e.g., at 9.9kms~^ at 1612 MHz, and 9.8 and 16.9 km s~^ 
at 1720 MHz (see Ta ble HI). Unfortunately, due the decay of the flux of XTE J1810-197 



( ICamilo et al.ll2007al ). we were not able to detect any OH absorption. 



Our OH absorption limits toward XTE J1810-197 are also meaningful for the OH seen 
at Vlsr ^ 28kms~^. The 1665 MHz OH feature at 30.7 km s~^ should have been detected at 
2-3(7 if it were between us and XTE J1810-197. However, we would not have expected to 
detect OH absorption at these velocities based on the H I absorption results (Table [2]). 



7. The Column Density to XTE J1810-197 



From the H I and OH observations that we have performed, it is possible to estimate 
the column density of hydrogen Nn, both atomic and molecular, to XTE J1810-197. This 
can then be compare d with the range of values i Vn = (6.5-14) x 10^^ cm~^ det ermined from 
fitting X-ray spectra (iGotthelf fc Halpernll2005l : iDurant fc van Kerkwijkl 120061 ). 



7.1. Estimate of the Atomic Column Density to XTE J1810-197 

We cannot directly measure the H I column density toward XTE J1810-197, since we 
do not have enough constraints to perform radiative transfer modeling in this direction. We 
can however still make an estimate of the column density to XTE J1810-197. To do this, 
we just integrate the H I emission spectrum between and 25.7 km s~^. Assuming that half 
of the emission at a given velocity comes from H I on the near side of the tangent point and 
the other half comes from gas with similar properties on the far side of the tangent point, 
we obtain an estimate of the column density to XTE J1810-197, 

1 r^s-^ 

AThi = - / 1.83 X 10^^ Tb{v) dv. (6) 

2 Jo 

Since there are clouds on the near side of the tangent point that significantly absorb the 
H I emission from the far side of the tangent point, as is evidenced by the absorption seen 
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against XTE J1810-197, this method provides a lower hmit to the column density of H I to 
XTE J1810-197. Upon performing the integration we find that A^hi ^ 1-8 ^ 10^^ cm~^. 



7.2. Estimate of the Molecular Column Density to XTE J1810-197 



We can use the pulsar "off" OH spectra to make an estimate of the molecular column 
density to XTE J1810-197. From Figures [3] and [6] we see that the 1612 MHz OH lines are 
in absorption while the 1720 MHz OH lines are in emission, with both having approximately 
the same amplitude. Such conjugate emission arises in regions where the OH column density 
is in the rang e 10^^ < Npyi/Av < 10^^ cni~^ kn i~^ s with Av the velocity resolution of the 
observations (jWeisberg et al.ll2005l : lElitzurlll992l ). The total column density of OH is found 
by integrating over the whole line. 

Integrating over the OH spectrum between and 25.7 km s~^ and assuming that half of 
the OH emission is from beyond the tangent point, we find 



(2.4 ± 0.3) X 10^^ < Nou < (2.4 ± 0.3) x 10^^ cm 



(7) 



Using standard abundances, the ratio of the numbe r of OH molecules to the number of 
hydrogen atoms is Nqh/Nh = 6 x 10~^ (lElitzurl 119921 ). which gives 



(4 ± 0.7) X 10^' < A^H < (40 ± 7) X 10^' cm 



(8) 



between us and the magnetar. 



Measur e ments of the ^"^00(2 -> 1) spectrum toward XTE J1810-197 are available from 
Dame et al.l ( l200lh . The total column density of hydrogen in molecular form can be deter- 
mined from the ^^C0(2 —>■ 1) spectrum using 

25.7 



ATh = 2iVH, = Xco / T^''{v)dv 



(9) 



where Xco is a conversion factor. The commonly used "standard" value for this is Xco = 
2 X 10^° cm~^ krn "^ s. However, Xco is known to vary depending on the line of sight 
( iMagnani et al.l Il998l ). The ^^C0(2 1) spectrum nearest to the line of sight toward 
XTE J1810-197 gives J^^'^ T§^{v)dv = 28.75 Kkms-\ of which we will assume that half 



arises from beyond the tangent point. From Figure 4 of iMagnani et al.l (1l998l ) we see 
that Xco = 10^° cm~^K~^km^^s is reasonable for the amount of ^^C0(2 1) toward 
XTE J1810-197. This then yields Nn > 2.8 x lO^^cm'^. 

This limit for obtained from the CO spectrum is consistent with the range of values 
determined using the OH spectra (eq. [8]). This also tells us that the OH column densities 
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are near the lower limit of Nq-h/Av = 10^^ cm ^ km ^ s. Since the CO spectra have a higher 
signal-to-noise ratio, we will use the CO value for the molecular gas contribution to A^h- 

Adding the molecular and atomic components of N-^, we obtain A^h ^ 4.6 x 10^^ cm^^. 
This limit is in agreement with the values determined from the X-ray spectrum of XTE J1810- 
197. Unfortunately we only have a lower limit for the total column density that is below the 
values determined from the X-ray spectrum. 



8. Models of the Line of Sight Toward XTE J1810-197 



A detail e d mo lecula r model of the ISM t o ward S OR 1806-20 has been developed by 
Corbel et al.l (119971 ) and ICorbel fc Eikenberryl (120041 ). depicted in Figure 8 of the latter. 
With increasing distance toward SGR 1806-20, their model encounters gas at Vlsr = 4, 
24, 30, 38, 44, and then 13kms ~^ in reaching the Scutum- Crux spiral arm (labeled as the 
30kms~^ spiral arm in Fig. 8 of ICorbel fc Eikenberryl |2004J ). Since the line of sight toward 
SGR 1806-20 is very close to that of XTE J1810-197 we might expect to encounter clouds 
at roughly the same velocities on the line of sight toward XTE J1810-197. 

The gas at 4km s~^ toward SGR 1806-20 is likely one of the two velocity components of 
the Heeschen Cloud ( iRiegel fc Crutcherlll972l ) and can be associated with the H I absorption 
feature at 7.7kms~^ toward XTE J1810-197, which is from the other velocity component 
of the Heeschen Cloud. The gas at 24kms~^ toward SGR 1806-20 can be associated with 
the H I absorption features at 22.8 or 25.7 km s"^ toward XTE J1810-197. We can associate 
the 13kms^^ gas toward SGR 1806-20 with the H I absorption seen at 14.1kms^^ toward 
XTE J1810-197. 



If the model of ICorbel et al.l (119971 ) and lCorbel fc Eikenberryl (120041 ) is correct, then we 
should also expect see H I absorption at velocities of roughly 30, 38, and 44kms~^ toward 
XTE J1810-197, since we observe absorption that can be associated with the 13kms~^ 
cloud in their model. In fact, we do not observe any H I absorpt ion against XTE J1 810- 
197 at these velocitie s . Thi s suggests that either (1) the model of ICorbel et al.l (119971 ) and 
Corbel fc Eikenberryl (120041 ) is not entirely correct, or (2) it cannot be applied to the line 
of sight toward XTE J1810-197, which is only 41' away, or (3) there is molecular material 
without H I along these lines of sight. The last possibility is not likely since molecular 
clouds are expected to have cosmic-ray ionization and photo-disso ciation regions with in and 
at their outer edges, which would produce atomic hydrogen (see iMinter et al.ll200ll . § 9.2, 
and references therein). 
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9. Discussion 



Using the DM = 178 ± 5 cm"^ pc measured fo r XTE J1810-197 flCamilo et all 120061 ). 
its distance according to tlie lCordes &: Lazid (120021 ) electron density model is 3.3 kpc. This 
model has a claimed average uncertainty of about 20% which, however, can be much larger 
for individual objects. For the sake of discussion, we assume an uncertainty of 1 kpc. The 
electron density model was derived using distances to pulsars that in many cases were de- 
termined via H I absorption spectra assuming a flat rotation curve, so that it seems most 
appropriate to compare the DM-derived distance of 3.3 ± 1 kpc with our fiat rotation curve's 
kinematic distance of 3.4to7kpc. These two v alues agree rema r kably well and imply that 



along the hne of sight to XTE J1810-197, the ICordes fc Lazid (120021 ) model gives a good 
representation of the average free electron density out to about 4 kpc. 

The distan ce to XTE J1810-197 has been estimated frora X-ray observations to range 
over 2.5-5 kpc (IGotthelf fc Halperru l2005l : iGotthelf et al.ll2004j ). These distances are deter- 
mined by converting A^h values obtained from fits to the X-ray spectra into visual extinction. 
Ay, and an estimate of the Ay per kpc in the Galaxy. This method is limited by a number of 
complications: (1) the spectral model used to fit the X-ray data, e.g., two blackbodies versus 
a blackbody and a power law; (2) the A^h versus Ay relationship determined locally (within 
~ Ikpc) but used for large distances (> Ikpc); (3) the large deviations from the fitted A^h 
versus Ay relationship for any particular line of sight; and (4) the d-Ay relationship also 
determined locally but used for large distances. 



Durant fc van KerkwiikI ((20061) determined ci = 3.1 ± 0.5 kpc toward XTE J1810-197 
ass uming an X-ray iVn = 14 x lO^^cm"^, whi c h is h igher than any of the values fitted 
bv IGotthelf et al. J2OO4J) or IGotthelf fc Halperd J2005h . We consider the d-Ay relation of 
Durant fc van Kerkwijk ( 20061 ) using red clump stars in the line of sight to XTE J1810-197 to 
be an improvement, although it is still necessary to apply an X-ray-fitted value of A^h to this 
relation and then convert it into a visual extinction Ay. Using A^h = 6.5 x 10^^ cm~^ deter- 



mined by IGotthelf fc Halperd (120051 ) , which is arguably a l ower limit, the extinction toward 



XTE J1810-197 becomes Ay ^ 3-4.5 mag using Figure 3 of IPredehl fc SchmittI Jl995h. This 
then y ields a distance estimate of 2-3.5 kpc based on Figure 7 of iDurant &: van Kerkwijk 



The kinematic distances determined from the H I absorption measurements presented in 
this paper rely only 01 1 the model of Galac t ic rot ation used to convert the me a sured velocity 
into a distance. The IWeiner fc SellwoodI (119991 ) and lEnglmaier fc Gerhard (120061 ) models 
both give consistent results (see Table [5]). We prefer the kinematic distance of 3.1-4.3 kpc 
determined in this paper over the DM- and X-ray-derived distances, because our conversion 
of measured velocity to distance is more direct and better constrained than through these 
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other methods. 



In Figure 7 of iDurant &: van Kerkwijkl (120061 ) we see that the extinction vs. distance 
rises steeply between 3.0 and 3.5 kpc and that for Ay < 13 the distance can be hmited to 
be less than 4 kpc. The line of sight toward XTE J1810-197' cont ains the Infrared Dark 
Cloud (IRDC) GlO.74-0.13 Jsimon et all bood : ICarev et all boooh . IRDCs are thought 
to be places where Giant Molecular Clouds are either just starting to form massive stars 
or on t he verge of beginning to form stars. The sudden increase in opacity vs. distance 
seen by iDurant fc van Kerkwijkl (120061 ) is likely associated with this IRDC. We note that 
the highest opacity reg i ons o f the IRDC cover less than half the area on the sky that 
Durant fc van Kerkwijkl (120061 ) used to determine the visual extinction vs. distance toward 
XTE J1810-197. Also, IRDCs can have extremely high opacities such that on ly the edge 
of the cloud is seen even in the infrared ( Simon et al. 20061: Minter et al. 2001 ). Since the 
opacities can be very large, it is plausible that iDurant fc van Kerkwijkl (120061 ) may have 
underestimated the actual extinction vs. distance toward XTE J1810-197. 



Jaffe et all (119821 ) associated IRDC GlO.74-0.13 with ^^CO (1 2) emission centered 
at 32kms~^, imply i ng d < 3.8 j]n'^kpc for a flat rotation cu rve and d < 4:.4:~^Q^kpc for the 
Weiner fc SellwoodI (119991 ) and lEnglmaier &: GerhardI (120061 ) models. Our OH spectra show 
a large spectral feature at the same velocities that we also associate with the IRDC. All 
observed H I absorption features lie at velocities lower than those associated with the IRDC. 
This puts XTE J1810-197 no further than the front edge of the IRDC. 



We consider the 4 kpc upper limit from Figure 7 of lDurant fc van Kerkwijkl (120061 ) to be 
a hard upper limit on the distance to XTE J1810-197. This further constrains slightly the 
distance estimate of 3.1-4.3 kpc obtained from H I absorption measurements (§ 4.3). Overall, 
we can thus summarize that the distance to XTE J1810-197 is 3.5 ± 0.5 kpc. Together with 
the measured proper mo tion of the AXP, this results in a transverse velocity corrected to the 
LSR of 212 ± 35kms~^ (IHelfand et al.ll2007l ). a perfectly ordinary velocity among pulsars. 



We would like to thank Jay Lockman for many useful discussions and T. Dame who 
graciously provided us with the ^^C0(2 — > 1) spectrum toward XTE J1810-197. The GET 
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Table 1. Observations of XTE J1810-197 at GBT 



Date 


Time 


(UTC) 


Species 


2006 Jun 6 


07:30 


- 11:00 


HI 


2006 Jul 22/23 


23:00 


- 07:30 


OH 


2006 Aug 31/Sep 1 


23:30 


- 05:00 


OH 


2006 Sep 2 


00:00 


- 05:00 


OH 


2006 Sep 4 


00:00 


- 05:00 


OH 


2006 Sep 24/25 


19:45 


- 03:40 


OH 


2006 Oct 21 


17:00 


- 19:09 


OH 


2006 Oct 21/22 


20:55 


- 01:45 


OH 
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Table 2. Gaussian fits to the H I absorption lines toward XTE J1810-197 



r 


^LSR 


FWHM 




(kms~^) 


(kms""*^) 


2.15 ±0.09 


7.73 ± 0.06 


3.3 ±0.2 


0.5 ±0.1 


14.1 ±0.2 


1.2 ±0.4 


1.53 ±0.08 


19.1 ±0.1 


3.8 ±0.4 


1.2 ±0.1 


22.8 ±0.2 


2.1 ±0.4 


0.94 ±0.1 


25.7 ±0.2 


2.4 ±0.5 



Note. — See Figs. [T] and [2] for data 
upon which these fits are based. 
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Table 3. OH absorption limits toward XTE J1810-197 



V 


a,- 


(MHz) 




1612 


0.09 


1665 


0.10 


1667 


0.10 


1720 


0.10 



1 a limits. 
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Table 4. Gaussian fits to the OH opacity in the pulsar "off" spectra toward 

XTE J1810-197 





r 


^LSR 


FWHM 


(MHz) 




(kms~^) 


(kms~^) 


1612 


0.026 ±0.002 


9.9 ±0.2 


4.4 ±0.6 


1612 


0.097 ±0.002 


29.54 ±0.05 


5.4 ±0.1 


1665 


0.059 ±0.009 


28.7 ±0.5 


7.9 ±0.4 


1665 


0.10 ±0.01 


30.69 ±0.08 


3.9 ±0.3 


1665 


0.020 ±0.002 


39.9 ±0.3 


5.5 ±0.8 


1667 


0.045 ±0.005 


28.1 ±0.6 


9.7 ±0.6 


1667 


0.086 ±0.007 


31.11 ±0.08 


4.3 ±0.3 


1667 


0.018 ±0.002 


42.0 ± 0.4 


6.3 ±0.9 


1720 


-0.023 ±0.002 


9.8 ±0.1 


3.1 ±0.3 


1720 


-0.020 ±0.002 


16.9 ±0.2 


4.8 ±0.5 


1720 


-0.021 ± 0.004 


26.2 ±0.2 


2.3 ±0.5 


1720 


-0.063 ±0.002 


29.5 ±0.1 


3.8 ±0.3 


1720 


0.014 ±0.002 


44.4 ± 0.3 


3.0 ±0.6 
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Table 5. Distance estimates for XTE J1810-197 



Measurement 


Method 


d 


Refs. 






(kpc) 





DM 


model fCordes & Lazio 2002) 

Av ^ 1.5 - 2.0 magkpc"^ -^d 


3.3 ± 1 


1 


X-ray: blackbody + power law 


- 5.0 


2 


X-ray: two blackbodies 


Av 1.5 - 2.0 magkpc"^ ->d 


- 2.5 


3 


X-ray: blackbody -I- power law 


A^H — ^ Ay —>■ red clump stars -^d 


3.1 ±0.5 


4 


X-ray: two blackbodies 


iVu — > Ay —>■ red clump stars —>d 


2-3.5 


5 


H I absorption 


Flat rotation curve fBurton 1988) 


3.4l|];^kpc 


5 


H I absorption 


Weiner & Sellwood ( 1999) model 


4.0;°jkpc 


5 


H I absorption 


Englmaier & Gerhard (2006) model 


3.7 ±0.6 kpc 


5 


H I absorption 


Brand & Blitz (1993) model 


2.4 ±0.5 


5 



Referenc es. — (1) [Camilo et al.l toO^ : (2) ICotthelf et all tooi) : (3) ICotthelf fc Habern 
|2005h : (4) burant fc van Kerkwiiir|2006l ): (5) this work. 



Note. — See § 4 for a discussion of the various rotation curve models that we fit to the H I data. 
The X-ray fits done by other authors are discussed in § 9. Overall, our best distance determination 
comes largely from the direct H I measurements on XTE J1810-197, and is 3.5 ± 0.5 kpc (see § 9). 
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Fig. 1. — Top: The H I emission spectrum (pulsar "off" spectrum) toward XTE J1810-197. 
Bottom: The H I absorption spectrum against XTE J1810-197. 
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Fig. 3.— Top: The OH 1612 MHz pulsar "off" spectrum toward XTE J1810-197. Bottom: 
The OH 1612 MHz absorption spectrum against XTE J1810-197. 
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Fig. 4.— Top: The OH 1665 MHz pulsar "off" spectrum toward XTE J1810-197. Bottom: 
The OH 1665 MHz absorption spectrum against XTE J1810-197. 
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Fig. 5.— Top: The OH 1667MHz pulsar "off" spectrum toward XTE J1810-197. Bottom: 
The OH 1667 MHz absorption spectrum against XTE J1810-197. 
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Fig. 6.— Top: The OH 1720 MHz pulsar "off" spectrum toward XTE J1810-197. Bottom: 
The OH 1720 MHz absorption spectrum against XTE J1810-197. 
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Fig. 7. — The flat Galactic rotation model {solid curve) of iBurtoru (119881 ). The x-axis is 
the distance from the Sun along the line of sight of XTE J1810-197. The y-axis is the 
radial velocity of the gas. The dotted lines indicate the deviation from Galactic rotation of 
±7kms~^. The arrow at ^sr = 25.7 km s~^ indicates the highest velocity H I absorption 
seen toward XTE J1810-197 (see Table[2]). The shaded regions indicate the allowed kinematic 
distances from the flat rotation curve (see discussion in §§ 4.2 and 4.3 for interpretation of 
these figures). 
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Fig. 9. — The same asFig. [3, except that it uses the Galactic rotation model of 
Englmaier fc GerhardI (120061 ). considering potential due to a bar and spiral arms. 
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Fig. 10. — The sam e as Fig. [71 except that it uses the Galactic rotation model of 



Brand fc Blita (119931 ). derived from observat i ons o f H II regions. The points and error 
bars are derived from Fig. 2 of iBrand fc Blitzl (119931 ). The solid line is a quadratic fit to the 
points. 
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Fig. 11. — Comparison of the GBT H I spectrum [black line) with the full-resolution SGPS 
H I spectrum {red line) toward XTE J1810-197. The SGPS spectrum convolved to the same 
resolution as the GBT is shown as the green line. The top panel shows the full range of 
emission, while the lower panel is zoomed in to show the difference in line widths for the 
Heeschen Cloud absorption around 8kms~^. 



